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ABSTFACT 

The  Plasma  Spectrophotoneter  constructed  at  the  Naval  Postgraduate 
School  in  1968  was  modified  by  the  addition  of  a  photomultiplier  cooling 
device  to  improve  time  resolution  and  increase  the  accuracy  of  the 
computed  values  of  the  relative  intensities  of  the  nitrogen  spectrum 
in  the  magnetically  confined  nitrogen  plasma  column. 

The  relative  intensities  determined  were  used  to  analyze  the  vari- 
ation in  local  electron  temperature  of  the  singly  and  doubly  ionized 
nitrogen  atoms  with  respect  to  magnetic  field.  The  local  electron  temp- 
erature was  found  to  vary  with  increasing  magnetic  field, reaching  a  peak 
temperature  of  8.1  electron  volts  at  4200  gauss. 

Evaluation  of  the  relative  number  densities  of  the  two  ionic 
species  indicated  that  the  local  electron  temperature  of  the  column 
was  the  local  electron  temperature  of  the  singly  ionized  nitrogen  atom. 
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I.   INTRODUCTION 

This  investigation  is  a  continuation  of  the  analyses  of  the  char- 
acteristics of  the  highly  ionized  plasma  column  produced  at  the  Naval 
Postgraduate  School  Plasma  Facility  which  have  been  undertaken  since 
its  fabrication  in  1962.  The  principle  objective  of  this  study  was 
the  analysis  of  shock  waves  in  the  magnetically  confined  plasma  column 
by  means  of  the  dual-Ebert  type  spectrophotometer  constructed  at  the 
Plasma  Facility  in  1968  by  Orlicki  [1] .  The  spectrophotometer  was  to 
be  modified  to  improve  time  and  spatial  resolution  and  then  used  to 
determine  the  variation  of  electron  temperatures  and  relative  number 
densities  of  key  ionic  species  with  changes  in  the  machine  parameters, 
such  as  magnetic  field. 

The  key  ionic  species  of  the  unperturbed  plasma  were  identified  as 
the  singly  and  doubly  ionized  atom  (N+  and  N++)  and  the  singly  ionized 
molecule  (N__)  [1] .  These  species  are  the  basis  of  the  continued  study. 
The  identified  lines  are  listed  in  Table  I. 


Table  I.  Identified  Lines  of  the  Nitrogen  Plasma 


N+  Wavelength  (A) 
4124.100 
4133.654 
4145.759 
4601.490 
4607.166 
4613.887 
4621.405 
4643.106 


N-H-  Wavelength  (A) 
4097.309 
4103.370 
4200.000 
4514.900 
4518.200 
4535.100 


A  complete  reanalysis  of  the  plasma  in  the  unperturbed  state  had 
to  be  accomplished  first  because  of  several  modifications  to  the  Plasma 
Facility,  including  the  addition  of  a  ten  inch  diffusion  pump  to  the 
cathode  end  of  the  system.  Thus  the  local  electron  temperature  and 
the  relative  number  densities  of  the  two  known  nitrogen  atom  species 
and  their  variance  with  respect  to  changes  in  the  magnetic  field  were 
to  be  determined  to  identify  changes  brought  about  through  incorporation 
of  the  modifications. 
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II.   BACKGROUND 

A.   SPECTROSCOPIC  DIAGNOSTICS 

The  nitrogen  plasma  produced  at  the  Naval  Postgraduate  School 
Plasma  Facility  emits  sufficient  light  for  spectroscopic  observation 
and  is  also  sufficiently  dense  and  long-lived  so  that  the  assumption 
can  be  made  that  the  velocity  of  the  electrons  is  very  nearly  Max- 
wellian  at  any  instant  of  time  and  at  any  point  in  space.  The  local 
electron  temperature  can  then  be  treated  as  a  function  of  time.  Since 
electron  impacts  dominate  almost  all  collisional  rate  processes,  the 
local  electron  temperature  is  essentially  the  quantity  determined  in 
spectroscopic  temperature  measurements. 

The  electron  density  is  usually  involved  in  the  relations  used  to 
calculate  the  electron  temperature  from  observed  intensities,  and 
therefore  temperature  and  density  measurements  are  not  really  inde- 
pendent. An  attempt  to  find  a  self-consistent  solution  of  the  various 
relationships  between  the  directly  measured  spectroscopic  quantities 
and  densities  and  temperatures  was  made  in  this  investigation. 

Since  the  assumption  that  the  plasma  is  in  local  thermal  equilibrium 
(LTE)  has  been  made,  the  relationship  based  on  the  fact  that  densities 
in  various  excited  states  are  proportional  to  the  products  of  the  stat- 
istical weights  with  the  exponentials  of  the  negative  ratios  of  excit- 
ation energy  and  the  thermal  energy  kT  may  be  used  [2] . 

Quantitatively , 

kT  =  d  E'  "  E  3 

Jin (IX  g'f'/I'X'  gf) 

where  the  quantities  are  defined  as  follows: 

11 


1.  First  Reference  Line 

I  =  relative  intensity 

A  =  wavelength  in  meters 

g  =  statistical  weight  of  the  lower  state  (2J  +  1) 

f  =  absorption  oscillator  strength 

E  =  excitation  potential  in  joules 

2.  Second  Reference  Line  —  same  as  above  except  designated  by  prime. 
Excitation  potentials  and  wavelengths  may  be  obtained  from  the  multi- 

plet  tables  of  Moore  (1959)  [3] ,  or  Striganov  and  Sventitski  (1968)  [4]  . 
Precise  wavelengths  and  statistical  weights  can  be  found  in  almost  all 
cases,  but  this  is  not  true  of  absorption  oscillator  strengths.  Except 
for  one-electron  systems,  the  oscillator  strengths  given  in  the  litera- 
ture are  approximate,  with  estimated  errors  ranging  from  10  percent  in 
the  case  of  neutral  and  other  relatively  simple  spectra  to  uncertainties 
of  factors  of  two  or  three  for  more  complex  spectra  [2,5] . 

Relative  intensities  can  be  measured  with  an  accuracy  of  better  than 
10  percent,  unless  the  wavelength  separation  of  the  lines  is  vary  large. 
Errors  due  to  uncertainties  in  the  oscillator  strengths  are  usually  the 
greatest  source  of  error,  and  uncertainties  of  30  percent  in  the  argu- 
ment of  the  logarithm  are  not  uncommon.  The  uncertainty  above  is 
reflected  in  a  similar  temperature  error,  and  since  kT  is  usually  of  the 
order  of  the  largest  separation  between  excitation  energies  of  non  reso- 
nance lines,  one  can  hardly  improve  accuracy  by  using  resonance  lines. 
The  gain  in  separation  between  excitation  energies  is  generally  offset 
by  larger  errors  in  oscillator  strengths,  in  relative  intensity  measure- 
ments, and  with  difficulties  with  self  absorption  [2,6]. 

Applying  the  concept  of  relative  intensity  to  different  ionic  species 
to  the  equation  for  absolute  intensity  of  a  spectral  line  given  by 


12 


Robinson  and  Lenn  [6]  ,  and  using  the  electron  temperature  determined 
from  the  equation  on  page  11,  one  can  arrive  at  the  relative  density 
equation: 


*1  W2    lh\ 

N2  "        I2glf      \X2  ) 


-(E2   -  EX)AT 
e 


where  the  quantities  are  as  defined  earlier,  except  that  they  are  now 
for  two  different  species.  Uncertainties  from  the  only  approximate 
validity  of  these  equilibrium  relations,  from  experimental  errors,  and 
also  from  errors  in  oscillator  strengths,  as  mentioned  above,  are  freq- 
uently of  similar  magnitude  which  implies  that  even  with  utmost  care  it 
is  often  not  possible  to  reduce  errors  in  relative  densities  below  30 
percent . 

It  must  be  emphasized  that  several  important  criteria  are  applicable 
when  using  the  above  equations.  These  expressions  are  only  applicable 
when  the  arc  pressure  is  held  constant  [2] ,  and  certain  ratios  of  para- 
meters in  both  equations  must  be  limited  if  meaningful  results  are  to 
be  obtained  [2,6,7] .  The  intensity  ratio  must  be  on  the  order  of  1.0 
to  6.0,  the  difference  in  excitation  potentials  must  be  greater  than  3.0 
volts,  and  the  wavelength  separation  of  the  lines  must  be  small. 

It  has  been  assumed  that  the  electron  density  is  sufficient  to 
maintain  a  Maxwellian  Distribution  of  quantum  state  population,  and  if 
this  is  not  so,  and  radiative  depopulation  of  lower  energy  states  takes 
place  to  any  significant  degree,  then  the  population  of  states  having 
energies  above  kT  would  be  favored  and  the  excitation  temperature  indi- 
cated by  the  line  intensities  would  be  higher  than  the  actual  electron 
temperatures  [2]  . 
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B.   THE  NITROGEN  ATOM 

The  nitrogen  atom  has  seven  identified  ionization  states.  However, 
for  this  investigation,  only  the  singly  and  doubly  ionized  states  are 
of  concern.  Their  ionization  potentials  are  29.593  and  47.426  electron 
volts  respectively  [8] .  The  lines  which  were  chosen  to  study  are  listed 
in  Table  II  along  with  their  relative  intensities  and  transitions.  All 
of  the  transitions  listed  are  allowed  [3,4] .  Four  lines  of  the  singly 
ionized  and  three  lines  of  the  doubly  ionized  atoms,  all  marked  by 
asterisk  (*) ,  are  listed  among  the  more  intense  lines  known  in  the 
nitrogen  spectrum  [4] . 


Table  II.  Relative  Intensities  and  Transitions  of 
Applicable  Nitrogen  Lines 


N+ 


NH 


o 

Wavelength  (A) 

I 

Transition 

4124.100 

4 

3s  P  -  3p  S 

4133.654 

5 

ii 

4145.759 

6 

ii 

*4601.490 

11 

_  3^0  _  3^ 
3s  P  -  3p  P 

*4607.166 

10 

ii 

4613.887 

9 

n 

*4621.405 

10 

it 

*4643.106 

11 

ii 

*4097.310 

10 

3s  S  -  3p  P 

*4103.370 

9 

ii 

4200.00 

6 

0  2^0  ^  2^ 
3s  P  -  3p  D 

*4514.900 

7 

_  4D0  _.  4^ 
3s  P  -  3p  D 

4518.200 

3 

ii 

4535.100 

2 

3p4S  -  3d4P° 

The  state  of  ionization  in  the  plasma  is  due  to  many  factors.  The 
primary  factor  is  the  cathode  current  interaction  which  supplies  energy 
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to  thermionically  emit  electrons  which  collide  with  the  gas  column. 
The  higher  the  cathode  current,  the  greater  number  of  electrons  emitted 
and  consequently,  the  greater  the  number  of  collisions  the  greater  the 
energy  exchange.  The  electrons  removed  in  the  ionization  process  may 
recollide  with  the  cathode  and  cause  secondary  emissions  which  are 
capable  of  causing  further  ionizations.  The  cathode-anode  accelerating 
potential  makes  further  energy  available  to  the  thermionically  emitted 
electrons,  and  so  it  continues.  Other  ionizations  may  occur  within  the 
plasma  due  to  atomic,  molecular  or  electron  interactions  which  would  be 
functions  of  other  machine  parameters  such  as  magnetic  field. 

C.   THE  PLASMA  FACILITY 

The  experiment  is  carried  out  in  the  steady  state  plasma  column, 

13   -3 
with  plasma  density  of  the  order  of  10  '  cm  ,  provided  by  a  hollow 

cathode  reflex  arc.  This  plasma  is  confined  in  a  ten-foot  long  column 

along  the  axis  of  a  six  coil  solenoid  which  produces  a  magnetic  field  of 

up  to  10,000  gauss  and  is  homogeneous  to  within  2.5%  along  the  axis  of 

the  plasma  column. 

The  plasma  tube  is  a  ten  foot  long  series  of  four  inch  diameter 
pyrex  sections  with  two  access  ports  available  on  each  side.  The  "theta" 
pinch  section  is  located  approximately  thirty  inches  from  the  floating 
anode,  and  the  access  port  used  for  the  collection  of  spectroscopic  data 
is  62  inches  from  the  cathode  (Figures  1  &  2) . 

The  magnetic  field  is  provided  by  six  main  magnet  coils,  a  cathode 
magnet,  a  first  anode  magnet,  and  a  mirror  magnet  at  the  floating  end. 
All  magnets  are  mounted  coaxially  to  the  tube  and  all  are  individually 
variable . 


15 


I 


Figure  1 
Plasma  Facility  -  End  View,  Cathode 


Figure  2 
Plasma  Facility  -  Side  View 
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The  nitrogen  gas  enters  the  tube  at  the  cathode  and  flows  through 
several  centimeters  of  the  cathode  tube  of  tantalum  where  it  is  ionized 
by  the  high  current  discharge  which  runs  into  the  interior  of  the  tube. 
This  highly  ionized  gas  then  leaves  the  cathode  and  forms  the  beam. 
The  cathode  is  self  heated  and  thus  electrons  are  thermionically  emit- 
ted from  its  inside  walls  or  by  field  emission  across  a  very  thin  sheath. 
Electrons  which  collide  with  gas  molecules  before  returning  to  the  walls 
are  trapped  and  drawn  out  of  the  cathode  by  the  axial  potential  gradient 
of  the  arc  as  well  as  by  the  pressure  gradient  of  the  plasma.  The  col- 
lisions of  the  electrons  with  the  gas  molecules  produce  positive  ions 
and  electrons,  radiation,  and  possibly  seme  metastable  excited  atoms. 
Nearly  all  of  these  products  are  trapped  in  the  tube  and  return  to  the 
wall.  This  heats  the  wall  to  temperatures  around  3000  C,  and  results 
in  thermionic  emission  of  electrons  [9] .  The  cathode  is  heated  locally 
to  high  temperatures  and  when  cathode  current  is  held  to  60  amperes, 
the  cathode  life  is  around  ten  hours.  The  duration  of  the  plasma  column 
in  a  steady  state  was  from  a  minimum  of  a  few  minutes  to  as  long  as  four 
hours. 

Argon  is  used  to  establish  the  initial  plasma  column,  and,  when  it 
has  stabilized,  the  system  is  switched  over  to  nitrogen  and  the  argon  is 
turned  off.  The  transition  is  readily  observable  in  the  column  as  it 
smoothly  turns  from  blue  to  pink.  Data  is  taken  when  the  nitrogen  beam 
is  established. 

D.   PLASMA  SPECTRDPHOTOMETBR 

The  plasma  spectrophotometer  is  a  dual-ebert  type  photoelectric 
spectrometer.  It  was  designed  to  accept  a  beam  of  light  through  a 
single  input  slit,  set  at  300  microns,  split  that  beam  in  half  by  means 
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Figure  3 
Internal  View  of  Plasma  Spectrophotometer 
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of  a  beam  splitter,  and  pass  the  resultant  beams  to  plane  gratings 
which  diffract  the  beams,  after  which  they  are  focussed  by  spherical 
mirrors  into  two  independent  chambers  [1] .   (See  Figure  3.)  The  dif- 
fracted beams  exit  through  separate  slits  and  are  picked  up  by  an 
Eldorado  Differential  Photometer  Unit  and  recorded.  The  Eldorado  is 
used  because  it  contains  all  the  apparatus  necessary  to  perform  this 
experiment,  and  is  highly  accurate,  giving  a  high  degree  of  spectro- 
scopic resolution.  It  provides  the  photomultiplier  heads,  power  supply 
and  recording  apparatus,  and  is  also  capable  of  mounting  both  end 
window  and  side  window  photomultiplier  tubes. 

The  1P21  photomultiplier  tube  was  used  for  this  experiment  as  it 

was  the  most  suitable  for  the  nitrogen  lines  under  investigation.  It 

o 
is  a  side  window  tube  with  maximum  sensitivity  at  4000  A  and  with  over- 

o 
all  sensitivity  from  3100-5400  A.  At  1024  volts  supply  voltage,  this 

tube  gives  a  current  amplification  of  2  X  10  with  an  equivalent  noise 

-12 
level  of  only  6  X  10    lumens.  Two  other  tubes  are  available,  the 

o 
6903,  head-on  tube,  which  has  a  maximum  sensitivity  at  4000  A  and  over- 

o 
all  sensitivity  from  2000-5600  A,  and  the  7102,  also  a  head-on  tube, 

o 
with  a  maximum  response  sensitivity  at  8000  A  and  overall  sensitivity 

o  ,  5 

between  5600  and  9600  A.  The  6903  has  a  current  amplification  of  4  X  10 

-12 

and  an  equivalent  noise  level  of  6.7  X  10    lumen.  It  is  better  suited 

for  work  in  the  ultraviolet  region.  The  7102  has  a  current  amplification 

5  -10 

of  4.6  X  10  and  an  equivalent  noise  level  of  7.5  X  10    lumens,  and 

is  better  suited  for  work  in  the  near  infra  red  region  [10] .  Using  any 

or  all  of  the  aforementioned  tubes  gives  the  capability  of  covering  the 

o 
entire  spectrum  of  the  spectrometer,  2000-9600  A. 
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Setting  the  required  wavelengths  is  accomplished  through  the 

use  of  digital  counters  capable  of  locating  each  individual  setting. 

o 
There  is  1/10  A  between  successive  settings.  The  gratings  are  mounted 

on  an  expanding  lever-arm  screw  type  mechanism  which  turns  the  grating 

and  the  digital  counter  at  the  same  time.  Settings  can  be  made  and 

accurately  reproduced  based  on  a  prepared  calibration  curve  (Figure  4) . 

This  curve  was  rechecked  with  a  mercury  arc  lamp,  after  all  modifications 

were  applied,  and  no  changes  were  found. 
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PLASMA  SPECTROPHOTOMETER 
WAVELENGTH  CALIBRATION  CURVE  (1) 


(Number  of  Turns 
vs.  Wavelength) 


#  of  turns 
300t 


250 


200 


L50 


100 


(10000) 
5000 


(12000) 
6000 


Figure  4 
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III.   EXPERIMENTAL  PROCEDURES 

A.   MODIFICATIONS  TO  SPECTROPHOTOMETER 

Since  time  resolution  was  desired,  it  was  necessary  to  modify  the 
differential  photometer  unit  through  the  addition  of  a  photomultiplier 
cooling  system.  The  cooling  of  the  photomultiplier  tubes  was  necessary 
to  reduce  the  detrimental  effects  of  dark  current  and  external  noise 
[11,12] .  A  cooling  device  was  constructed  which  reduced  the  dark  cur- 
rent and  external  noise  pickup  rapidly  and  conveniently  through  the 
use  of  liquid  nitrogen.  No  difficulty  was  experienced  with  electrical 
short  circuits,  but  the  usual  problem  of  window  fogging  had  to  be 
resolved . 

A  schematic  drawing  of  the  cooling  device  is  seen  in  Figure  5.  The 
side  window  photomultiplier  tube  housing  from  the  Eldorado  unit  was  wrap- 
ped with  a  coil  of  copper  tubing  ending  in  the  exit  lens  housing.  The 
original  aluminum  lens  housing  was  replaced  with  a  nonconducting  material, 
polyvinylchloride  tubing  (PVC) ,  which  had  the  proper  internal  dimensions 
and  was  machinable.  A  gas  vent  tube  was  installed  in  the  top  of  the  PVC 
tubing  to  prevent  pressure  buildup  inside  the  system,  and  to  allow  water 
vapor  to  escape.  The  entire  copper  coil  was  connected  to  the  tube  housing 
with  Thermon,  Heat  Transfer  Cement  [13] .  Since  the  tube  housing  was 
aluminum,  welding  or  brazing  was  impractical,  and  the  heat  transfer  cement 
insured  effective  heat  conduction  between  the  metals.  The  entire  system 
was  then  covered  with  polystyrene  for  further  insulation. 

In  operation,  liquid  nitrogen  is  placed  in  a  reservoir  connected  by 
a  Tee  device  to  the  input  end  of  both  coils.  The  liquid  nitrogen  fills 
the  coils  to  the  height  of  the  liquid  in  the  reservoir  and  then  begins 
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evaporating.  Immediate  cooling  of  the  tube  housing  takes  place  through 
conduction.  Because  the  walls  of  the  housing  cool  first,  water  vapor 
that  is  present  condenses  upon  them,  thus  preventing  harmful  conden- 
sation on  the  tube  or  its  wiring.  The  nitrogen  gas  flows  upward  through 
the  coil  and  fills  the  lens  housing,  cooling  the  photomultiplier  tube  and 
the  lens.  Excess  gas  escapes  through  the  vent  tube.  To  prevent  fogging, 
a  1/16  inch  coil  of  nichrome  wire,  six  inches  long,  and  wrapped  in  a 
fiberglass  sheath  to  prevent  radiation,  was  wrapped  around  the  circumfer- 
ence of  each  lens,  and  connected  to  a  variac.  Adjustment  of  the  variac 
as  needed  produces  sufficient  heating  of  the  lens  to  prevent  fogging. 

Thermocouples  inserted  in  the  base  of  each  photomultiplier  tube  hous- 
ing register  the  temperature  of  the  resistor  bank  of  each  tube.  The  photo- 
multiplier tubes  are  cooled  sufficiently  for  use  fifteen  or  twenty  minutes 
after  filling  with  liquid  nitrogen.  About  two  liters  of  liquid  nitrogen 
is  needed  to  cool  the  tubes  to  about  -80  C,  and  about  1/4  liter  of  liquid 
nitrogen  should  be  added  each  10  minutes  for  a  thirty  minute  period.  Data 
can  be  taken  for  about  30  minutes  after  this  without  further  cooling.  The 

dark  current,  with  1012  volts  across  the  tube,  is  reduced  by  a  factor  of 

3 

about  10  within  20  minutes  of  applying  liquid  nitrogen  to  the  reservoir. 

There  is  essentially  no  difference  in  the  spectral  sensitivity  of  the  tube 
at  low  temperatures  [10] . 

B.   OPERATION  OF  THE  SPECTROPHOTOMETER 

The  spectrophotometer  is  rolled  into  position  opposite  the  port  to 
be  used  in  the  Plasma  Facility.  The  unit  is  raised  and  leveled  so  that 
the  input  optics  are  positioned  as  required  [1] .  The  power  switch  of 
the  Eldorado  Unit  is  activated  and  one  liter  of  liquid  nitrogen  is  placed 
in  the  cooling  reservoir.  The  variac  is  turned  on  and  adjusted  to  heat 
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the  output  lenses.  After  fifteen  minutes  have  elapsed,  another  liter 
of  liquid  nitrogen  is  added.  It  takes  30  minutes  for  the  electronic 
equipment  to  settle  out  and  after  that  interval  the  high  voltage 
switch  is  activated.  A  lens  cover  is  placed  over  the  input  optics 
and  the  meters  for  both  photomultipliers  are  calibrated  [14] .  After 
removal  of  the  lens  cover,  data  may  be  taken. 

The  first  line  of  interest  is  selected  by  use  of  the  digital  count- 
ers and  the  input  optics  are  adjusted  to  give  maximum  intensity.  The 
"B"  channel  photomultiplier  intensity  reading  is  used  as  the  base 
intensity  in  all  cases,  and  the  "A"  channel  intensity  reading  is 
adjusted  to  equal  "B".  The  method  known  as  "leap-frogging"  was  used 
exclusively  in  this  investigation ,  i.e.,  starting  with  the  same  line 
on  each  channel,  with  equal  intensity  on  both,  and  moving  on  to  the 
next,  registering  equal  intensities  on  both  channels  before  moving  on 
to  the  next  [1] .  Thus  all  intensities  end  up  being  referenced  to  the 
first  line  to  come  under  investigation.  In  this  investigation,  it  was 
always  the  longest  wavelength  under  study. 

C.   EXPERIMENTAL  DATA  AND  COMPUTATIONS 

A  nitrogen  plasma  column  is  established  at  the  Plasma  Facility  and 
the  parameters  of  the  system,  magnetic  field,  arc  current,  arc  pressure, 
and  anode  voltage  are  recorded.  The  magnetic  field  desired  is  selected 
and  held  constant  while  all  the  selected  nitrogen  lines  are  investigated, 
It  takes  approximately  30  minutes  to  evaluate  14  lines.  The  magnetic 
field  is  then  changed  and  the  parameters  are  checked  for  deviation. 
This  is  continued  until  all  the  lines  have  been  evaluated  at  every 
magnetic  field  in  question.  For  this  investigation,  the  magnetic  field 
was  set  at  2400,  3000,  3600,  4200,  and  4800  gauss.  The  arc  current  was 
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constant  at  60  amps  and  the  anode  voltage  varied  from  64  to  75  volts 

during  separate  runs.  The  arc  pressure  varied  during  any  particular 

-4  -5 

run  from  1.6  X  10   to  7.6  X  10   torr.  This  will  account  for  seme 

variance  in  the  results. 

The  intensities  of  each  line  are  recorded  as  shown  in  Figure  6, 
and  the  uncorrected  relative  intensities  (URI)  are  computed.  The  URI 
of  line  1  is  the  reading  for  line  1,  as  both  channels  have  been  balanced 
with  channel  "B"  being  the  base.  The  URI  for  the  remaining  lines  is 
computed  as  follows:  The  URI  of  line  2  is  to  the  URI  of  line  1  as 
the  average  reading  for  line  2  is  to  the  reading  of  channel  "B"  (base 
channel)  for  line  1.  Any  changes  in  readings  for  the  same  line  on  tavo 
channels  (caused  by  moving  the  gratings)  must  be  taken  into  account, 
and  thus  the  use  of  the  average  reading  for  line  2  above.  The  compu- 
tations are  continued  substituting  line  2  for  line  1  and  line  3  for 
line  2  proceeding  through  the  remainder  of  the  lines  in  the  same  manner. 
The  URI  must  then  be  corrected  for  spectral  sensitivity  of  the  photo- 
multiplier  tubes  used.   (See  Figure  7) .  The  URI  multiplied  by  the 
correction  factor  gives  the  Corrected  Relative  Intensity  [1] . 

A  computer  program  was  written  to  calculate  the  temperature  equation 
given  in  the  "Spectroscopic  Diagnostics"  portion  of  this  paper.  How- 
ever, before  temperatures  could  be  computed,  the  values  of  the  maximum 
and  minimum  intensity  ratios,  the  maximum  wavelength  separation,  and 
the  minimum  excitation  potential  difference  had  to  be  determined.  The 
program  was  adjusted  so  that  one  parameter  could  be  varied  while  holding 
the  others  constant,  and  was  broad  enough  so  that  the  parameter  varied 
could  assume  all  values  in  question  during  one  run  of  the  computer. 
Limits  on  the  values  of  parameters  were  established  so  that  consistant 
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PLASMA  SPECTROPHOTOMETER  SAMPLE  DATA 


CHANNEL  "A" 

CHANNEL  "B" 

WAVE- 

Ionic 

APPROX. 
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APPROX. 
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1P2 1  SPECTRAL  SENSITIVITY  AND  RESULTANT  CORRECTION  FACTORS 


1P21  SPECTRAL  SENSITIVITY  CURVE 
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values  of  temperature  were  obtained.  Three  different  sets  of  relative 
intensities  determined  on  three  separate  runs  of  the  Plasma  Facility 
were  used  to  keep  the  diagnosis  unbiased.  Limits  established  for  the 
parameters  were  maximum  relative  intensity  of  4.0,  minimum  relative 
intensity  ratio  of  0.25,  maximum  wavelength  separation  of  500  angstroms, 
and  minimum  excitation  potential  difference  of  7.0  volts.  The  computer 
program  with  all  corrections  included  and  all  parametric  varying  devices 
removed  is  shown  on  page  36. 

A  computer  program  to  calculate  relative  number  density  ratios  of 
the  two  nitrogen  species  was  written  and  modified  with  the  inclusion 
of  the  above  determined  parameters  [1] ,  and  is  shown  on  page  39.  The 
temperature  used  in  this  computation  must  be  taken  from  the  temperature 
calculation  above  and  entered  as  a  constant.  The  average  temperature  is 
usually  chosen,  however  any  temperature  within  the  standard  deviation  of 
the  mean  will  suffice. 
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IV.   RESULTS  AND  CONCLUSIONS 

A.   RESULTS 

1.  Electron  Temperatures 

Entering  data  similar  to  that  shown  in  Figure  6  into  the 
electron  temperature  program  led  to  quantitative  results  for  the 
electron  temperature.  The  local  electron  temperature  was  found  to 
vary  with  increasing  magnetic  field  and  is  shown  in  Figure  8  for  the 
case  of  the  singly  ionized  nitrogen  atom. 

The  temperature  varied  from  3.85  electron  volts  at  2400  gauss 
to  a  peak  temperature  of  8.1  electron  volts  at  4200  gauss.  It  drop- 
ped rapidly  to  4.4  electron  volts  at  4800  gauss.  The  compression 
of  the  plasma  column  into  a  smaller  volume  as  the  magnetic  field  is 
increased  can  be  physically  observed  from  the  side  of  the  Plasma 
Facility.  The  density  of  the  singly  ionized  atoms  must  then  increase 
and  consequently  atomic  interactions  increase  giving  rise  to  higher 
temperatures.  Above  4200  gauss,  the  electron  temperature  decreases 
rapidly  plausibly  due  to  the  collision  rate  increasing  at  such  a  rate 
as  to  cause  a  corresponding  decrease  in  diffusion  rate,  and  thus  a 
greater  distribution  of  energy  throughout  the  column.  These  temp- 
eratures, Figure  8,  closely  approximate  those  found  by  Major  Thomas  J. 
Haycraft,  USA,  whose  studies  of  the  plasma  column  under  the  same  con- 
ditions were  conducted  using  probe  techniques  and  frequency  analyses. 

The  temperatures  determined  for  the  doubly  ionized  nitrogen 
atom  were  not  usable  as  they  were  consistently  unrealistic  and  ir- 
regular. Temperatures  varied  from  10  electron  volts  to  237  electron 
volts,  well  above  arcing  potential.  Evaluation  of  individual 


30 


ELECTRON  TEMPERATURE  VS  MAGNETIC  FIELD 
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temperatures  showed  that  out  of  six  temperatures  determined  for  each 
machine  run  at  constant  magnetic  field,  four  were  negative  temperatures. 
Investigation  of  the  lines  giving  these  data  showed  that  each  was  an 
allowed  transition,  and  that  two  of  the  lines  are  noted  as  among  the 
more  intense  nitrogen  lines.  Improper  identification  of  lines  is  not 
the  problem,  as  in  each  case,  there  are  no  nitrogen  lines  or  argon  lines 
which  could  be  mi  scons  trued  as  the  identified  lines  [4]  .  Their  relative 
intensities  are  not  consistent  with  those  relative  intensities  shown 
in  Table  II,  but  the  relative  intensities  of  the  singly  ionized  species 
are  consistent.  There  is  no  readily  acceptable  explanation  for  these 
deviations  because  the  relative  intensities  for  all  lines  are  determined 
in  a  single  run  with  a  single  base.  As  for  the  negative  temperatures, 
it  seems  to  follow  that  the  doubly  ionized  nitrogen  atom  is  not  in  local 
thermal  equilibrium,  and  that  the  velocity  of  the  electrons  is  not  Max- 
wellian.  The  electron  density  is  probably  such  that  radiative  depop- 
ulation of  lower  energy  states  is  prevalent,  giving  higher  electron 
temperatures  than  the  actual  temperature. 
2.  Relative  Number  Density,  N++/N+ 

It  was  necessary  following  the  results  of  the  electron  tempera- 
ture calculations  to  assume  that  the  local  electron  temperature  for  the 
singly  ionized  nitrogen  atom  was  the  local  electron  temperature  for  the 
system  for  the  purpose  of  obtaining  relative  number  densities  of  the 
two  species  studied.  The  average  local  electron  temperature  for  the 
singly  ionized  nitrogen  atom  for  a  particular  magnetic  field  was  entered 
as  a  constant  in  the  relative  number  density  program,  and  the  correspond- 
ing computational  data  were  entered  and  the  program  run.  The  data  and 
the  temperature  were  changed  for  each  value  of  magnetic  field  studied 


32 


33  A 


RELATIVE  NUMBER  DENSITY  VS 


MAGNETIC  FIELD 


29 


25 


21 


N-H- 
N+ 


17 


13 


5 


2.4 


3.0 


3.6 


4.2 


4.8 


B    (Kilogauss) 
Figure  9 


33 


(See  Figure  9) .  As  the  magnetic  field  increased  to  4200  gauss,  the 
relative  number  density  decreased  in  accordance  with  the  postulated 
increase  in  the  density  of  the  singly  ionized  atom.  The  density  of 
the  doubly  ionized  atom  must  also  decrease  somewhat,  as  the  temperatures 
calculated  show  an  approach  towards  the  local  equilibrium  state,  i.e., 
from  48  electron  volts  to  18  electron  volts  at  3600  gauss.  This  does 
not  in  any  way  assume  that  the  condition  of  local  thermal  equilibrium 
is  achieved  by  the  doubly  ionized  nitrogen  atom. 

B.   CONCLUSIONS 

This  investigation  has  identified  the  variance  in  local  electron 
temperature  with  magnetic  field  for  the  singly  ionized  nitrogen  atom 
and  determined  that  the  singly  ionized  nitrogen  atom  is  in  local  thermal 
equilibrium  in  the  plasma  column.  The  variance  in  temperature  with  mag- 
netic field  of  the  doubly  ionized  nitrogen  atom  was  not  determined  to 
a  satisfactory  degree.  The  doubly  ionized  nitrogen  atom  was  determined 
not  to  be  in  local  thermal  equilibrium. 

It  was  shown  that  the  local  electron  temperature  of  the  system  was 
the  local  electron  temperature  of  the  singly  ionized  nitrogen  atom  as 
evidenced  by  the  comparison  of  local  electron  temperatures  to  relative 
number  densities  determined  at  the  same  temperature  (Figure  10) . 

As  the  "theta  pinch"  device  became  operational  late  in  the  investig- 
ative period,  and  still  requires  further  modification,  it  was  not  possible 
to  include  an  analysis  of  shock  wave  perturbations  on  the  nitrogen  plasma 
column.  It  is  felt  that  this  can  be  accomplished  in  the  future  utilizing 
the  equipment,  as  modified,  and  the  procedures  developed  in  this  investi- 
gation . 
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ELECTRON  TEMPERATURE  COMPUTER  PROGRAM 


,2X,F7.2,2Xt«PERCENT« 

21  DC  80 ' 1=1 ♦ J 

THE  ABOVE  DO  LOCP  CHECKS  TO  SEE  IF  ALL  CATA  HAS  EEEN  ENTERED 


READ 

fcPITE  (6.12) 

IF  (J)  100,100,21 
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READ    (5,2)     EXCIFC( I ),JLPPER(I ) ,JLCW( I ) ,F(  I  ),WL(  I  ), 
1RELINK  I) 
WFITE(6tlll    WL(I) 

8C  CONTINUE 
THE  ABOVE  INSTRUCTION  WILL  READ  DATA  CN  A  SERIES  CF  "J" 
SPECTRAL  LINES. 

N  =  1.0 

WRITE  (6,3) 

DC  22    1=1, J 

GLCW(I)=(2.0*JLCW(I))*1.0 

WRITE  (6,4)  N,WL( I),EXCIPC( I ),JUPPER( I),JLCW( I), F( I), 
1RELINTU  )  vGLQUdl 

N=N+1.0 

22  CONTINUE 

DC  502  Z=l,3 

AFMAX(Z+1)=ARMAX( ZJ-1.0 

APMIN(7)=1.0/ARVAX(Z+1) 

R£TMAX=ARMAX(Z+1 ) 

RATMIN=ARFIN(Z) 

N  =  1.0 

S=l 

WRITE  (6,5) 

DC  2  6  1=1, J 

DC  66  MAP=I,J 

KPIS=EXCIPC(MAP)-EXCIFC(I  ) 

CHECK=ABS(KRIS) 

IF  (CHECK. LT. VE )GC  TC  70 
THE  ABCVE  INSTRUCTION  CHECKS  TO  SEE  IF  THE  DIFFERENCE  IN  EX- 
CITATION POTENTIAL  IS  GREATER  THAN  7.0  ELECTRCN  VOLTS.   THE 
DIFFERENCE  MUST  BE  SEVERAL  VCLTS  FCR  VALIC  RESULTS. 

WAVF=ABS(WL( I )-WL( VAP) ) 

IF  (WAVE.GT.VWL)GC  TC  70 
THE  ABOVE  INSTRUCTION  CHECKS  TO  SEE  IF  THE  DIFFERENCE  IN 
WAVELENGTHS  IS  LESS  THAN  500  ANGSTRC^S.   MUST  BE  FCR  VALID 
RESLLTS. 

IF  (KRIS)  24,  IOC, 23 

23  P£TIO(S)=RELINT(I)/RELINT(MAP) 
VANG=PATIC(S) 

IF  ( VANG.GT.RATVAXJGC  TC  70 

IF  ( VANG.LT.RATMMGC  TC  70 

TEMP(S)=(((EXCIFO(MAP)-EXCIPC(I) )*(EVCONV) )/B)/ALCG 
1( (PATIO(S) ) *(GLCW( MAP) /GLOW ( I ) ) * ( F ( MAP ) /F ( I ) )* ( WL ( I )** 
23.0/WL(MAP)**3.C) ) 

GO  TO  25 

24  RATIC(S)=RELlNT(MflP)/RELlNT(I) 
VUNG=RATIC(S) 

IF  ( VUNG.GT.RATVAXJGC  TC  70 
IF  ( VUNG.LT.RATVINJGC  TC  70 

TEMP(S)  =  (  ( (EXCIFO( I  )-EXCIPC(MflP) )*(EVCCNV ) )/B )/ALCG 
1( (PATIO(S) )*(GLCW( I )/GLCW(MAP) )*(F(I)/F(MAP)  )*(WL(MAP ) 
2**3.C/WL(I)**3.C) ) 
GC  TO  2  5 
7C  TEMP(S)=0.0 
R£TIC(S)=0.0 
THIS  PROGRAM  IS  SET  LF  TC  RECCPC  A  VALUE  CF  0.0  FCR  TEMPERA- 
TURE ANC  RATIO  IF  THE  EXCITATION  POTENTIAL  DIFFERENCE  OR 
WAVELENGTH  DIFFERENCE  DC  NCT  MEET  ABOVE  CRITERIA.   IT  WILL 
ALSO  ADJUST  THE  VALUE  CF  N  FOR  THE  CCMPUTATION  CF  THE  MEAN, 
ETC.. 

25  WFITE  (6,6)  N , TEM F ( S ) , WL ( I ) , WL ( V AP ) , R AT  10 ( S ) 
N  =  N-H.O 

S  =  S+1 
6  6  CONTINUE 

26  CONTINUE 
SLM=0.0 
ASLM  =  C0 
N  =  N-1 
S=S-1 

DC  27  I=1,S 

TEMPd  )=ABS(TEVF(  I  )  ) 

IF  (TEMP(I) ,E0.C.C)N=N-1.0 

27  SLM=SUV+TEMP(I ) 
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IF  (N.LE.l. 
AVF-=SUM/N 
IF  (N.EQ.l. 
DT  28  1=1, S 
IF(TEMP( I ). 
PIFF=ABS(AV 
SG=DTFF**2 
Gf    TC    128 

123    OIFF=0.0 
SC=0.0 

126  ASLM=ASUM+S 
SIGM=SGRT( A 
SIGP=SIGM/S 
WRITE  (6,20 
WFITE  (6,20 
WFITE  (6,20 
WRITE  (6,14 
WRITE  (6,10 
PAVE=AVE*SI 
QAVE=AVE-S! 
WFITE  (6,21 
PEPR=SIGP*1 
WRITE  (6,21 
EAVE=AVE*P/ 
ESIGP=SICB* 
WRITE  (6,20 
PEAVE=EAVE+ 
GEAVE=EAVE- 
WFITE    (6,21 

502    CCNTINUF 
GC    TO    20 

ICC    STCP 
END 


0)GC    TC    5C2 

0)GC    TC    5C2 

EQ.O.OJGC    TC    123 
E-TE*P(I  )  ) 


SUM/(N-1.0) ) 
G*T(M 

0)  VfcL 

1)  VE 

2)  RATVAX,RATMN 
)     N 

\     AVE,SIGR 

GB 

C  R 

0)  PAVE, CAVE 
OO.C/AVE 

2)  FERR 

EVCCNV 

8/EVCCNV 

5)  EAVE,ESIGP 

SIGP 

SIGR 

1)  FEAVE,CEAVE 


SAMPLE  CATA 


26 
28 
26 
21 
21 


36C 
361 
3*2 
07C 
C62 


21.C61 
21.C6C 
21.C*4 

3C. 331 
3C. 33C 
3S.640 
36.250 
3  fl  .  2  3  0 
41.51C 
-ICC 


J  ( U » 

2.0 
2.0 
2.0 
2.0 
1.0 
1.0 
0.0 
1.0 

1.5 

0.5 
2.5 
3.5 
0.5 
1.5 


CAPCS 

. 

. 

J(L) 

1. 

c 

2. 

0 

3. 

c 

1. 

0 

0. 

c 

1. 

0 

1. 

c 

2. 

0 

0. 

5 

0. 

5 

1. 

2. 

5 

0. 

5 

1. 

5 

(NCT 


C    SCALE): 


WAVELENGTH 


KREL) 


118 

+4124.1CCE- 

■1C 

21.49 

lie 

♦4133.654E- 

-1C 

35.07 

118 

♦4145.759E- 

-1C 

41.52 

133 

+46C1.490E- 

■1C 

96.65 

320 

♦4607.167E- 

■1C 

89.90 

0798 

♦4613.387E- 

■1C 

60.23 

106 

+4621.405E- 

-1C 

89.44 

0796 

♦4643.106E- 

■10 

107.50 

484 

♦4097.3CSE- 

-1C 

13.74 

Z<*2 

+4103.37CE- 

•10 

10.46 

398 

♦4200.000E- 

■1C 

157.05 

286 

♦4514.900F- 

■1C 

12.68 

178 

♦4518.200E- 

-1C 

52.06 

304 

♦4535. 100E- 

■1C 

22.04 
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RELATIVE  NUMBER  DENSITY  COMPUTER  PROGRAM 


THIS    PROGRAM    CALCULATES    THE    RELATIVE    NUMBER    DENSITIES    CF    TWO 
SPECIES    IN    A    PLASMA    FRO*    THEIR    SPECTRAL    LINE    DATA.       THE    TWO 
SPECIES    APE    CALLED    Nl    AND    N2,WITH    THE    N2    SPECIES    HAVING    THE 
HIGHER    EXCITATICN    POTENTIAL     (IN    GENERAL)     AND    BEING    ENTERED 
AS    THE    SECOND    OF    TWO    DATA    SETS. 

DIMENSION    GLOW (A, 50) ,  F(4,50) , EXC I PO (4 , 50 ) ,  PEL  INT ( 4,50 ) 
1WL(4,50) , JLOW(4,50) , JUPPER (4 , 50 ) ,  RAT  10 ( 50, 50  ) , 
2NUMDEN(50,50) 

REAL    GLOW,F,EXCIPC,RELlNT,B,WL,TEMP,EVCONV,RATIO,JLOW, 
1JUPPER,KRIS,N0NE,NTWC,NUMDEN,N,MATE 

EVCONV=1.30207E-1<9 

B=1.38042E-23 

VVL=500.E-10 

VE=7.0 
THE  APOVE  2  CARDS  GIVE  THE  MAXIMUM  WAVELENGTH  ANC  MINIMUM 
EXCITATION  POTENTIAL  DIFFERENCES  USED  IN  THE  COMPUTATIONS. 

RATMAX=^.0 

RATMIN=0.25 
THE  TWO  CARDS  ABOVE  SET  A  MAX  AND  A  MIN  CN  THE  RELATIVE 
INTENSITY  PATIO  OP  ANY  TWC  SPECTRAL  LINES. 

TFMP=44700.0 
TEMPERATURE  IS  CONSIDERED  A  CCNSTANT  FCR  THESE  CALCULATIONS. 
THE  VALUE  OF  TEMPERATURE  USED  SHOULD  BE  OBTAINED  FRCM  THE 
ELECTRON  TEMPERATURE  CCMFUTER  PROGRAM  FOR  THE  SPECTRAL  CATA 
USED  HFPEIN. 

1  FORMATt •l'.^X, 'ELECTRON  TEMPERATURE  =  ',F9.0,2X, 
l'DEGPEES  KELVIN*) 

2  FORMAT  (15) 

3  FORMAT  ( '1  •  ,20X,'CATA  CHECKS1/////) 

4  FORMAT  (4F10.0,E13.0,7X,F10.0) 

5  FORMAT  (20X,4PE13.3/) 

6  FORMAT  ( «1  •  ,3X, 'NO. »,5X, 'WAVELENGTH' ,7X, 'EXCITATION 
lPCTENTIAL'UlX,  'J',17X,'J',  11X,  • ABSORPT ION • ,4X. 

2 'RELATIVE ',5X,'2J+1'/13X,' (METERS)' , 1 5X, •  ( VCLTS) •  ,  1 IX  , 
3MUPPER  LEVEL)  ',5X,  MLGUER  LEVEL )', 6X ,' STR ENGTH ', 5X, 
4'  INTENSITY' ,4X  ,  •  ( LCW )  • /// ) 

7  FORMAT  (  1X,F4.0,F4.0,1X,4PE13.3,12X,OPF7.3, 16X.F4.1, 
114X,F4.1,10X,F8.5,6X,F7.2,5X,F4.1/) 

8  FOPMAT  ( '1  ',4X,'N2.',3X, • WAVEL ENGTH ( 2 ) • ,5X  ,  • Nl . • , 3X, 
1  'WAVELENGTH( 1)  ',10X,  '  INTENSITY  RATIO', 11X, 'NUMBER 
2DENSITY  RATIOS' /13X,' (METERS) ', 16X, • (METERS) • ,15X, 
3'(I(2)/K1))',1<3X,'N(2)/N(1)'///) 

9  FORMAT  (/4X,F4.0,3X,4PE13.3,4X,0PF4.0,3X,4PE13.3, 12X, 
11PE10.3,20X,2PE10.2) 

10  FORMAT  ( '1' ,14X,« SUM  =   •  ,4PE12 .2,// 15X, • AVER  AGE 
1NUMBER  DENSITY  RATIC  =   •  ,2PE10 . 2//// / 15X, • A8S (NBAR-N ) 
2'  ,15X, »ABS( (NBAP-N)SO) • ) 

11  FORMAT  (/15X,2PE10.2,18X,2PE10.2) 

12  FORMAT  ( *1  •  ,14X, 'THE  TCTAL  NUMBER  OF  ACCEPTABLE 
1CCMPUTATICNS  (I.E., WITHIN  SPECIFIED  LIMITS)  MACE  =   ', 
2F5.0) 

13  F0RMAK///15X,  'STANCARD  DEVIATICN  =   •  ,  12X  ,  E  1C.4/  /15X 
1, 'STANDARD  DEVIATICN  CF  THE  *EAN  =   »,E10.4) 

14  FORMAT  (///////////20X,'THE  NUMBER  DENSITY  RATIO  OF 
1SPECIES  2  TO  SPECIES  1  =   • ,2PE 10.2 , 2X ,'♦/-', IX, 0PE8. 2 

200  FOPMAT  ( /////10X,  'MAXIMUM  WAVELENGTH  CIFFERENCE  =   ', 
13FE10.1,2X, 'METERS' ) 

201  FORMAT  (/////10X, 'MINIMUM  EXCITATICN  PCTENTIAL 
1DIFFFFERENCE  =   •  ,F4.  1 ,  2X ,' VCLTS •///// ) 

202  FORMAT  ( 10X ,  'MAX  I  MUM  ALLOWABLE  INTENSITY  RATIC  =  ', 
1F6.1/////10X, 'MINIMUM  ALLOWABLE  INTENSITY  RATIC  =  ', 
2F5.2///I 

SPECIES  2  =  N++,  AND  SPECIES  1  =  N*  FCR  THIS  INVESTIGATION. 
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5CC 

5C2 


1 

THE  AB 
TCTAL 

THIS  P 
30 


1 

THE  AB 
TCTAL 

THIS  F 
31 


1 

FCR  TH 
THE  0A 

32 


33 


THIS  I 
CREATE 
VALID 


THIS  I 
OR  LES 


1 
2 

34 

IF  THE 

ENCE  D 
DENSIT 
35 
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READ  (5.2) 
IF  (INDIA) 
WPITE  (6,3) 
DP  30  1=1,1 
READ  (5,4) 
WL(1,I ).REl 
OVE  CARD  WI 
OF  ••INDIA" 
WPITE  (6,5) 
RINT-OUT  CH 
CONTINUE 
READ  (5,2) 
WRITE  (6,3) 
DC  31  1=1, J 
READ  (5,4) 
WL(2,I  ),PFL 
OVE  CARD  WI 
OF  "JULIET" 
WPITE  (6,5) 
RINT-OUT  CH 
CONTINUE 
NCNF=1.0 
NTW0=1.0 
WRITF  (6,6) 
KILO=INDIA 
DC  33  K=l,2 
DC  32  1=1, K 
GLCW(K,I )=( 
WRITE  (6,7) 
JLOW(K,I) ,F 
E  RESULTING 
TA  SET  AND 
NTWO=NTWO-H 
CCNTINUE 
NCNE=N0NE+1 
NTW0=1.0 
KIL0=JULIFT 
CCNTINUE 
N0NE=1.0 
NTW0=1.0 
WRITE  (6,8) 
DC  39  K=1,J 
DC  38  1=1,1 
KPIS=EXCIPr 
CHECK=ABS(K 
IF  (CHECK. L 
NSURES  THAT 
R  THAN  OP  E 
RESULTS. 
MATE=ABS( WL 
IF  (MATE.GT 
NSUPES  THAT 
S  THAN  500 
PATIC(K,I )= 
VCNG=PATIO( 
IF  (VCNG.GT 
IF  (VONG.LT 
NUMDEN(K.I  ) 
F(2,K)  )*( (W 
(B*TEMP) )  )  ) 
GC  TO  35 
NUMDEN(K.I  ) 
PATIO(K,{)= 
EXCITATION 
0  NOT  MEET 
Y  RATIO  ANT 
WFITE  (6,9) 
NUVPEN(K,I ) 
NCNF=N0NE+1 
XPAY^NONE 
CONTINUE 
NTW0=NTWC+1 


I  NOI  A 

501  ,501,502 

NDIA 

EXCIFO(l.I)  tJUPPERUtl  )  ,JLCW(  1,  I),F(1,  I), 

INT(l.I) 

LL  REAC  THE  M  CATA  SET  CONSISTING  OF  A 
NUMBER  CF  CARCS. 

WL(  1,1  ) 
ECKS  TC  SEE  THAT  ALL  CATA  IS  ENTEREC. 

JULIET 

ULIET 

EXCIP0(2,I) ,JUPPER(2,I ),JL0W(2,I),F(2,I), 

INT(2,I) 

LL  REAC  THE  N2  DATA  SET  CONSISTING  OF  A 

NUMBER  OF  CARDS. 

WL(2,I) 
ECKS  TC  SEE  THAT  ALL  DATA  IS  ENTEREC. 


ILC 
2.0 

NO 
(K, 

PR 
THE 
.0 

.0 


*JLCW(K,I ) )+1.0 

NE,NTWC,WL(K,I ),EXCIPC(K,I),JUPPER(K,I), 
I) ,RELINT(K,I)  ,GLCW(K,I) 

INT-CUT,  THE  FIRST  NUMBER  WILL  INDICATE 
SECCND  WHICH  ONE  OF  THAT  DATA  SET. 


ULIET 

NDIA 

(2,K)-EXCIFC(1,I) 

RIS) 

T.VE)GC  TC  34 

THE  DIFFERENCE  IN  EXCITATION  PCTENTIALS 
QUAL  TC  7.0  VOLTS  WHICH  IS  NECESSARY  FCR 


IS 


(2,K 
.VWL 

THE 
ANGS 
RELI 
K,I  ) 
.RAT 
.RAT 
=  (RA 
L(2, 


)-WL(l,I ) ) 

DIFFERENCE  IN  WAVELENGTHS  IS  EQUAL 
TPGMS.  NECESSARY  FCR  VALIC  RESULTS 
NT(2,K)/RELINT(1,I  ) 


TO 


VAX)GC    TC    34 

flMGC    TC    34 

TIC(K,I)*(GL0W(1,I )/GL0W(2,K)  )*(F(  1.  I)/ 

K)/WL(1,I))**3.0)*(EXP( (KRIS*EVCONV)/ 


=  0.0 
0.0 

POTENTIAL  CIFFERENCE  CR  WAVELENGTH  DIFFER- 
ABOVE  CRITERIA,  THE  VALUES  OF  THE  NUMBER 

THE  INTENSITY  RATIO  WILL  eE  SET  TC  ZERO. 
NTWC,WL(2,K) ,N0NE,WL(1, I ) , R AT  I C (K ,  I  ) , 

.0 


.0 
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3<5 


ZF 

NC 

CC 

N  = 

SU 

AS 

DC 

DO 

NL 

IF 

THE  ABOV 

TIONS  ON 

CCRRECT. 

45  SU 

46  CC 
AV 
WR 
DC 
DC 
IF 
DI 
SO 
GC 
DI 
SO 
WR 
AS 
CC 
SI 
SI 
WF 
WR 
WR 
WR 
WR 
WR 
WR 

PRIN 
THE 


BPA=NTWC 

NE=1.0 

NTINUE 

(XRAY-1.0)*( ZEBRA-l.O) 

M=0.0 

UM=0.0 

46    K  = 

45  1  = 
MDEMK 

(  NUWD 
E  CHEC 
LY,     IN 


1,  JULIET 
I, INDIA 

,1  )=ABS(NUFDEMK,I)  ) 
EMK.I  ).EC.O.O)N=N-1.0 

K  CORRECTS  THE  COUNT  TO,  ACCEPTABLE  CCMPUTA- 
CRDER  TO  KEEP. THE  STATISTICAL  DISTRIBUTION 


I  ) 


AVE 


TO  123 


123 

12A 
47 

48 


THE 
♦  /- 


5C1 


M=SUM4-NLMQEN(K 

NTINUE 

E=SUM/N 

ITE  (6,10)  SUV 

48  K=l , JULIET 

47, 1=1, INDIA 

(NUMDSN(K,I  ).EC.O.0)GC 
FF=ABS(AVE-NUVDEN(K,I)) 
=DIFF*2 

TC  124 
FF=0.0 
=  0,0 

ITE  (6,11)  DIFF,SC 
UM=ASUM+SQ 
NTINUE 

GM=SQFT(ASUM/(N-1.0)  ) 
GB=SIGM/SQRT(M 
ITE  (6,12)  N 

(6,13)     SIGM,SIGB 
(6,1)     TE^F 

(6.200)  VWL 

(6.201)  VE 

(6.202)  RATI*AX,PATMIN 
(6,1*1     AVE,SIGB 

T-OUT  ABOVE  GIVES  THE  AVERAGE  NUMBER 
STANDARD  DEVIATICN  CF  THE  MEAN. 
OP 
D 
.  • 


SAMPLE  DATA  CARCS  (NCT  TC  SCALE) 


ITE 
ITE 
ITE 
ITE 
ITE 
ITE 


DENSITY  (N2/N1) 


ST 
EN 


2S 

26 
2B 


E 

8 
36C 
361 
,362 
21.070 
21.062 
21.061 
21.060 
21.06* 

6 
3C.331 
3C.330 
3S.64C 
38.250 
38.230 
41. 510 
-ICO 


J  fill  J(L)        F  WAVELENGTH  I(REL) 

2.0  1.0  .118  4-4124. 100E-1C  21.49 

2.0  2.0  .118  +4133.654E-10  35.07 

2.0  3.C  .118  +4145.759E-1C  41.52 

2.0  1.0  .133  4-4601. 490E-1C  96.65 

1.0  0.0  .320  +4607.167E-1C  89.90 

1.0  1.0  .0798  +4613. 887E-1C  60.23 

0.0  l.C  .106  +4621.405E-1C  89.44 

1.0  2.0  .0796  +4643.106E-10  107.50 

1.5  0.5  .484  -1-4097. 309E-1C  13.74 

0.5  0.5  .242  4-4103. 370E-10  10.46 

2.5  1.5  .398  +4200.000E-10  157.05 

3.5  2.5  .286  +4514.900E-1C  12.68 

0.5  0.5  .178  4-4518. 200E-1C  52.06 

1.5  1.5  .304  4-4535. 100E-IC  22.04 
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